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We study inclusive heavy quark and exclusive heavy quark pair produc- 
tion in pp, pA and AA interactions. Intrinsic transverse momentum is 
introduced in pp interactions. Nuclear effects, limited to kj- broadening 
and nuclear shadowing, are introduced in pA and AA interactions. The 
nuclear dependence is studied over a range of energies, both in fixed 
target and collider setups. 



Introduction 

Previously, we calculated the heavy quark production cross section in pp interac- 
tions. In that work, we focussed on the scale dependence and fixed the quark mass 
and scale parameters by comparison to the total cross section data [||]. In this exten- 
sion to pA and AA interactions, we study the effects of shadowing and kx broadening 
on the spectra of heavy quark production. 

* This work was supported in part by the Director, Office of Energy Research, Division of Nuclear 
Physics of the Office of High Energy and Nuclear Physics of the U. S. Department of Energy under 
Contract Number DE-AC03-76SF00098. 



We first briefly update the calculation of the total cross section with new parton 
distribution functions and then go on to discuss the relevant differential distributions. 



Total cross section 



Recall that at leading order (LO) heavy quarks are produced by gg fusion and qq 
annihilation while at next-to-leading order (NLO), qg + qg scattering is also included. 
At any order, the partonic cross section may be expressed in terms of dimensionless 
scaling functions f-j'''^ that depend only on the variable r] [Q] 

a,,{s,ml,f,') = [Ana^if^yfj: fl^^{^ In' f4) > (1) 

where s is the partonic center of mass, mg is the heavy quark mass, yU is the scale, 
and 7] = s/ArrLQ — 1. The cross section is calculated as an expansion in powers of 
as with A; = corresponding to the Born cross section at order 0{a1). The first 
correction, k = 1, corresponds to the NLO cross section at 0{a^). It is only at this 
order and above that the renormalization scale dependence enters the calculation since 
when k = 1 and / = 1, the logarithm ln(/i^/mQ) appears. The dependence on the 
factorization scale, the argument of as, appears already at LO. We assume that the 
renormalization and factorization scales are the same. The next-to-next-to-leading 
order (NNLO) corrections to next-to-next-to-leading logarithm have been calculated 
near threshold p[ but the complete calculation only exists to NLO. 
The total partonic cross section can also be written as 



crij{s,mQ, fi"^) = J didus^ 



d aij{s, t, u) 



dtdu 

where the double differential partonic cross section at leading order is proportional 
to the Born cross section, 

,d^c^f^{s,i,u) 
didu 



''Zl ' =Sis-i-u- 2ml)af^{s, t, u) . (3) 



The expressions for the Born cross section can be found in e.g. Ref. [Q. The NLO 
and NNLO cross sections in the qq channel are also proportional to a^, albeit with 
more complicated scaling functions. In the gg channel, the higher order cross sections 
are again proportional to the Born cross section to leading logarithm and as are the 
highest powers of the factorization scale dependent logarithms 0. However, beyond 
leading logarithm, the scaling functions also depend on the color state in which the 
QQ pair was produced. The expressions for the Born cross section in both channels 
may be found in Ref. along with the NLO and NNLO expressions near threshold. 

The cross sections may be calculated in either one-particle inclusive kinematics or 
pair invariant mass kinematics, depending on whether one or both heavy quarks are 
detected. Single heavy quark hadroproduction in pp interactions is defined by 

p(Pi) + p(P2) ^ Q{pi) + X{px) , (4) 



where X denotes any hadronic final state containing the heavy antiquark and Q{pi) 
is the identified heavy quark. The reaction in Eq. (D is dominated by the partonic 
reactions 

qih) + qih) QiPi)+Xmp2), (5) 

g{k,)+g{k2) Q(pi)+X[g](p2). (6) 

At LO or if (P2) = ^(^2)? the reaction is at partonic threshold with Q momentum 
P2. At threshold the heavy quarks do not have to be produced at rest but are instead 
produced with equal and opposite momentum, leading to trivial values of the pair 
transverse momentum {pt^ +Pr2 ~ 0) azimuthal angle (0 = vr). The Mandelstam 
invariants on the hadronic level are 

s = iP, + P,Y , t=iP^-p^r , u = {P,-p,)\ (7) 

with the corresponding partonic invariants 

S={h + k2f , i={k2~pif , u={ki-p^f. (8) 

The pair invariant mass kinematics are somewhat different. Now the pair is treated 
as a single unit and we have 

p{Pi)+p{P2) ^ QQ{P')+X{Px). (9) 

On the parton level, Eqs. and are now 

q{kr) + q{k,) QQ{p') + X{px) ^ (10) 

g{k,)+g{k2) QQ{p') + X{px) . (11) 

The square of the heavy quark pair mass is p'"^ = M^. In the case when px 7^ 0, 
more invariants than the three typical Mandelstam invariants can be defined 0]. 
At partonic threshold, px = 0, only the three Mandelstam invariants remain with 
s = and 

t-ml = - — {l-(3,,cosd), (12) 
u-rriQ = — — (1 + /9m cos6') , (13) 



where (3m = yl ~4mQ/M2, and 9 is the scattering angle in the parton center of 
mass frame. 

The total hadronic cross section is obtained by convoluting the total partonic cross 
section with the parton distribution functions of the initial hadrons. 



app[s,mQ} 



/ -t)F[{xi, fi^)F'f'{x2, fi'^)aij{T,mQ,fi^) , (U) 

■ - „ J irrii^/s T 



where the sum i is over all massless partons and Xi and X2 are the hadron momen- 
tum fractions carried by the interacting partons. The parton distribution functions. 



denoted by Ff , are evaluated at the factorization scale, assumed to be equal to the 
renormalization scale in all our calculations. 

In our previous work on pp interactions we obtained our results with two sets of 
parton distribution functions that had successfully predicted the general behavior of 
the low X deep-inelastic scattering data P], MRS D— ' Q and GRV HO 0. Much more 
low X data has become available since then which shows that these parton densities 
are too high at low x and new fits to the most recent data [§] are now available |0. 
(The low X region is particularly relevant for QQ production at the LHC.) Our latest 
calculations employ the most recent NLO parton distribution functions from the MRS 
group, MRST HO (central gluon) |TD|. This group has also made a leading order set 
available We compare the MRS D-', MRST HO, and MRST LO results in Fig. |l| 
and in Tables |T] and 0, all updated from our earlier work |]T| . Table |I| has also been 
expanded to include the fixed target energies at which we compare the pp and pA 
calculations. As before, all our numerical results are obtained using the NLO code 
by Mangano, Nason and Ridolfi (MNR) §. 

All the cross sections are calculated at LO and NLO with the same mass and scale 
parameters obtained from the pp study — nic = 1.2 GeV, fic = 2mc, rub = 4.75 GeV, 
and /ife = rrih. The agreement with the total charm cross section data is at least as 
good with the MRST sets as with the MRS D— ' set. In fact, the somewhat larger 
NLO cross section at fixed target energies with MRST HO tends to agree better with 
the data than the older set, in part because some of the gluon density at low x has 
been redistributed to higher x. However, at the pp and AA collider energies shown in 
Tables |l| and |^, the growth of the cross sections is considerably slower with the new 
distribution functions than with the MRS D— ' set. 

The MRST LO results are for comparative purposes only. A LO set is not con- 
sistent with a NLO calculation just as a NLO set is inconsistent with a strictly LO 
calculation. We note however that the tabulated results are all calculated with as 
at two loops, regardless of whether the LO or NLO result is given, because in the 
MNR code is calculated to two loops regardless of the order of the chosen parton 
density. Since the QCD scale. A, is tuned to fit all available data in global analyses of 
parton distributions, whether at LO or NLO, the number of loops in the evaluation 
is important. We remark however, that it is rather standard procedure to evaluate 
all orders of the hadronic cross section with the same NLO parton density and a 
two loop calculation of a^. This is done so that one can study the variation of the 
partonic cross section alone while keeping the nonperturbative results fixed. Also, 
beyond NLO, no standard NNLO or resummed evaluations of the parton densities 
exist, generally because the partonic matrix elements are limited to near threshold 
beyond NLO. In addition, the data used in the global analyses of the parton densi- 
ties are typically far from threshold and of only limited use in such an evaluation. 
Therefore the uncertainties in such an analysis would be large. 

We note that when a fully LO calculation of the charm cross section is done with 
a one loop evaluation of as in the global analyses of the parton densities, as with the 
MRST LO set, the LO cross sections in Table |l] increase by ~ 60%. This change 
is almost solely due to the difference between the one and two loop evaluations of 
as- The MRST LO A is 0.204 GeV when Uf = 3, leading to a^-^""? = 0.287 and 



^2-ioops ^ Q for /i = 2mc. The A associated with the MRST HO set is larger, 
A = 0.353 GeV when nj = 3, corresponding to aJ-^°°P = 0.364 and 0"^'^°°^' = 0.263 
at the same scale. Thus a fully LO calculation with MRST LO partons would be in 
better agreement with the tabulated LO results for MRST HO partons. These results 
show that one must use caution when using a LO set with a NLO calculation or a 
NLO set for a LO calculation. 
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Figure 1: Total charm production cross sections from pp and pA measurements ||T^, 
|T^, |I^, |T3|, |TB| compared to calculations with nic = 1.2 GeV and /i = 2mc- The 



parton distributions used are MRS D-' (solid); MRST HO (dashed) and MRST LO 
(dot-dashed). Updated from Ref. 



Differential cross sections 
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Table 1: Total cc production cross sections at fixed-target and collider energies. 
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Table 2: Total hh production cross sections at collider energies. 



While the total cross section predicts the yield of heavy quark production over all 
phase space, it cannot provide much useful information on nuclear effects such a pt- 
broadening and shadowing. It is also independent of the choice of kinematic schemes. 
Any kx broadening will not affect the total cross section but will have a strong 
influence on the pt distributions. Shadowing may reduce or enhance the nuclear 
cross section relative to that of the proton but the effect may be more apparent in 
some regions of phase space than others. To obtain more information on nuclear 
effects, it is thus necessary to turn to distributions. In addition, a real detector does 
not cover all phase space. The differential distributions can be tuned to a detector 
acceptance. 

The double differential hadronic cross sections are written in the factorized form 

2 (i (Tpp{s,t,u) s--^ dxi dX2 j.p, 2\£Pf 2\ n- I- X -\-2'^ ^«i('^' ^' ^) /i r-N 
S = 1^1 / fi{xi,fi )ff{x2,^i )jK{s,t,u)s ^j— 15 

dtdu ■ „ Jx- xi Jx- X2 dtdu 

where /f is the parton density, /f(a;, /x^) = F[{x, fi'^)/x, and J'k is a kinematics- 
dependent Jacobian. These results are difficult to calculate analytically and expres- 
sions only exist for s'^d'^aij{s,t,u)/dtdu, see e.g. 0. However, the MNR code cal- 
culates single and double differential distributions as well as total cross sections. In 
our first paper |l|, we neglected any fragmentation effects on the produced heavy 
quark and intrinsic transverse momentum of the initial partons because the inclusive 
D meson pt distributions agreed with the bare charm quark distributions, showing 



no need for such effects. However, as pointed out in Ref. iT^, the pair distributions 



particularly the pair pt and the azimuthal angle distributions at fixed target energies 
are broader than can be accounted for by the bare quark distributions. The general 
behavior of the pair data was reproduced in Ref. [O] by including a Peterson-type 



fragmentation function |T8[ along with an intrinsic kx of the initial partons. Indeed, 
the combination of fragmentation with an intrinsic average of 1 GeV^ leads to 
single D meson px distributions in good agreement with the bare charm distribution. 
These effects can be added by modifying Eq. (|1^) to include the additional integrals 



/ dz,dz,d^kr,d^kT. ^^/^^''' ^^/Q^"^' ^'\ ,{kTMkT2) . (16) 

J Z4 



The fragmentation functions Dijiq{z) describe the hadronization of the heavy 
quark Q into hadron H. The final-state hadron H carries a fraction z of the heavy 
quark momentum and is assumed to be coUinear to the parent heavy quark. The 



Peterson fragmentation function [jT8|, fit to D production data in e^e collisions, 

Dh/q{z) = — —, 77- TT^ , (17) 

z{l-l/z-eQ/{l-z)y 



is used in the code. The parameter eg is 0.06 for charm and 0.006 for bottom |]T9 . 
The normalization N is defined by J2h I DH/Q{z)dz = 1. 

The Peterson function results in a nearly 30% decrease in the average charm quark 
momentum during hadronization. This decrease must be compensated for in inclusive 
production because the decreased momentum of the D meson would no longer agree 
with the single D px distributions. The intrinsic fc-r of the initial partons is such a 
compensation mechanism which has been used successfully to describe the pt distri- 
butions of other processes such as Drell-Yan pO[ and hard photon production pl| . 



Some of the low px Drell-Yan data has subsequently been described by resummation 
to all orders but the inclusion of higher orders has not eliminated the need for this 
intrinsic kx- It was found that assuming {k"^) = 1 GeV^ in pp and 7cp interactions 
produces agreement with the inclusive data on the level of that found with the bare 
quarks alone and does a better job of reconciling the exclusive DD data with the 
calculated distributions |T7 . 



The implementation of intrinsic kx in the MNR code is not the same as in other 
processes due to the nature of their calculation. The cancellation of divergences 
is handled numerically in the code. Since adding additional numerical integrations 
would slow this process, the kx kick is added in the final, rather than the initial state. 
In Eq. (p!6D, the Gaussian function gp{kx), 

dvi^T) = exp(-4/(4)p) , (18) 

Ti\Kxlp 



with (A;|-)p = 1 GeV^ fl^, multiplies the parton distribution functions, assuming the 



X and kx dependencies completely factorize. If this is true, it does not matter whether 
the kx dependence appears in the initial or final state modulo the caveat described 



below. In the code, the QQ system is boosted to rest from its longitudinal center 
of mass frame. Intrinsic transverse momenta of the incoming partons, kxi and /ct2, 
are chosen at random with fcf,^ and distributed according to Eq. (jl^). A second 
transverse boost out of the pair rest frame changes the initial transverse momentum of 
the QQ pair, p-^ in Eqs. ([10|) and (|ll|), to + kxi + /cr2- The initial fcy of the partons 
could have alternatively been given to the entire final-state system, as is essentially 
done if applied in the initial state, instead of to the QQ pair. There is no difference 
if the calculation is to LO only but at NLO an additional light parton also appears 
in the final state. In Ref. |2^, it is claimed that the difference in the two methods is 
rather small if < 2 GeV^. The effect of the intrinsic fcy decreases as the center of 
mass energy increases. 



broadening in nuclei 

The average intrinsic kx is expected to increase in pA interactions. This broaden- 
ing is observed in Drell-Yan [^, J/ip ||2^, and T |]23| production and has been used 



to explain high pion production in nuclear interactions p5| . Since kx broadening 
has not been explicitly measured in charm production, we follow the formulation of 
Ref. ||25|, similar to Ref. ||26[ for J/ip and Drell-Yan production where the kx broad- 



ening arises from multiple scattering of the projectile partons in the target, 

(4)A = (4)p + ((^)-l)A'(/i) • (19) 

We define the number of collisions in a proton-nucleus interaction, averaged over 
impact parameter, as 

. , Jd%Tl{b) 3 , , 

where TA^b) = J dzpA{b,z) is the nuclear profile function, and use this to calculate 
the average broadening. In the definition of (z/), a^N is the inelastic nucleon-nucleon 
scattering cross section, po is the central nuclear density, and Ra is the nuclear radius. 
Our calculations are done with A = 200, Ra = 1.2^4^/^, and po = 0.16/fm^. 

The strength of the nuclear broadening, A^, depends on p, the scale of the inter- 
action 

Thus is larger for bb production than cc production. This empirically reflects the 
larger kr broadening of the T |2^ relative to the J/ip |^. We evaluate A^(yu) at 
H = 2mQ for both charm and bottom production. We find ((z/) — l)A^(/i) = 0.35 
GeV^ for charm and 1.57 GeV^ for bottom in pA collisions at 6 = and A = 200. We 
can change the centrality by changing (u) — 1. Reducing (z/) — 1 by two corresponds 
to a more peripheral pA collision while multiplying (z/) — 1 by two corresponds to 



central AA collisions. We will use the central values given above in our pA results 
and the larger values, 0.70 GeV^ for charm and 3.14 GeV^ for bottom, in our AA 
results. Note that the total {k^)A in both cases is obtained by adding {k^)p = 1 GeV^ 
to the above nuclear broadening. 

The effects of broadening alone on several pair and single inclusive quantities 
are shown for a fixed target beam energy of 158 GeV in Figs. ^ and ^ The changes 
in the shape of the pair quantities in Fig. |^ are rather striking. The bare charm 
pair distributions alone are compared to results with (k^) = 1, 1.175, 1.35, and 1.7 
GeV^. The differences in the pair pt and distributions are quite apparent. The 
change in the pair pt distribution is due entirely to kr broadening since fragmentation 
effects would steepen the solid curve. The broadening is still important at large pair 
Pt, causing more than an order of magnitude difference in the cross section between 
(k^) = 1 GeV^ and 1.7 GeV^ at = 4 GeV. Likewise, introducing intrinsic kx 
changes the distribution from sharply peaked at = vr for the bare charm quarks 
to nearly flat when {k^) = 1 GeV^ to a peak at = when (A;|,) = 1.7 GeV^. On 
the contrary, the pair rapidity shows no effect of /ct broadening. 

The broadening effect on the single charm distributions shown in Fig. ^ is not 
as strong. Note that all the distributions, including the charm rapidity distribution, 
are affected by broadening and fragmentation. The charm pt distributions for the 
bare quark and including fragmentation with (k^) = 1 GeV^ are almost identical, in 
accord with the observation that no intrinsic kr was needed to describe the D meson 
Pt distributions. The difference between the bare quark and (k^) = 1.7 GeV^ is less 
than a factor of five for pt = 3.5 GeV. Fragmentation is the most important effect 
on the xp and rapidity distributions since there is little change in the distributions 
as (k^) is increased. Any type of fragmentation, even that of a delta function, will 
cause the rapidity and xp distributions to be modified. This is because fragmentation 



enhances these distributions in the central region ||2^. If Dh/q{z) is a delta function 



the momentum does not change, pn = Pq, but Ej^ = Eq — nxq + mjj, resulting in a 



rapidity shift |^ with 



dpzQ dp^H mT,H coshyndyH coshy ndyn 

dn oc dyg = — = — — = ^ —i^^^=^^^= (22) 

Eq Jrri^^H cosh^ yn - mjj + mj^ J cosh^ yn - 



where 



2 _ ml -ml 



a' = " ^ ^ . (23) 
mrpfj 

For mc = 1.2 GeV, m^) = 1.87 GeV and mp^D ~ V^rriD, cP' = 0.25, increasing the 
D cross section at = by ^ 15%. The total cross section is unchanged. When 
m[, = 4.75 GeV, ms = 5.27 GeV and mT,B = V^^b, «^ = 0.09, increasing the 
B cross section by ^ 5% at ys = 0. The range of the enhancement is {ynl < 2.5, 
independent of energy. If the Peterson function is used instead, increases to 






Figure 2: Effects of kx broadening alone on exclusive NLO cc pair production in 
pA interactions at 158 GeV as a function of pt, y, and 4>. The per nucleon cross 
section is given. The curves are the bare quark distributions alone {{k^)A — and 
no fragmentation) (sohd) and, including fragmentation, {k'^) — 1 (dashed), 1.175 
(dot-dashed), 1.35 (dotted) and 1.7 (dot-dot-dot-dashed) GeV^. 



increasing the D cross section at T/i? = by ~ 30% for {z) ~ 0.7 and the B cross 
section by ~ 15% at ?/s = for {z) ^ 0.85. These (z) values are typical for the 
Peterson function with the eg values given above. The xp distribution is also affected 
because xp = 2mr sinhy/i/i. 
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Figure 3: Effects of fc-r broadening alone on inclusive NLO c quark production in 
pA interactions at 158 GeV as a function of pt, xp, and y. The per nucleon cross 
section is given. The curves are the bare quark distributions alone {{kp)A = and 
no fragmentation) (solid) and, including fragmentation, (k^) = 1 (dashed), 1.175 
(dot-dashed), 1.35 (dotted) and 1.7 (dot-dot-dot-dashed) GeV^. 



Nuclear shadowing 



Experiments p8[ have shown that the proton and neutron structure functions 
are modified in a nuclear environment. For momentum fractions x < 0.1 and 0.3 < 



X < 0.7, a depletion is observed in the nuclear parton distributions. The low x, or 
shadowing, region and the larger x, or EMC region, is bridged by an enhancement 
known as antishadowing for 0.1 < x < 0.3. For simplicity, the entire characteristic 
modification as a function of x has been referred to as shadowing. 

Shadowing is an area of intense study with numerous models available in the 
literature However, none of the models can satisfactorily explain the behavior 
of the nuclear parton distributions over the entire x and range. Therefore, we 
choose to use parameterizations of shadowing based on nuclear DIS data. As in DIS 
with protons, the nuclear gluon distribution is not directly measured and can only be 
inferred from conservation rules. We use the recent EKS98 parameterization based on 
the GRV LO parton densities in the proton. (See the contribution to this volume 
by K.J. Eskola et al. for further details.) Each parton type is evolved separately above 



yUo = 1-5 GeV |2^, 0. The initial gluon distribution shows significant antishadowing 
for 0.1 < X < 0.3 while the sea quark distributions are shadowed. Figure ^ compares 
shows the shadowing ratios at ^ = and /i = 10 GeV. The valence ratio changes 
little with scale but the evolution of the gluon and sea ratios with scale is more 
significant. The gluon ratio changes somewhat more rapidly with scale than the sea 
ratio but the difference is not large. Less shadowing is seen in the antishadowing 
region at larger scales for the sea quark ratio while the gluon ratio decreases with 
scale in the same region. The effects of shadowing alone on the DD and BB decay 
contributions to the dilepton continuum have been studied in Ref. . 



Shadowing should depend on the spatial location of the parton in the nucleus. 
Unfortunately, there is little data on this topic. One relevant experiment, Fermilab 
E745, studied the spatial distribution of nuclear parton densities with uN interac- 
tions in emulsion. The presence of one or more dark tracks from slow protons is used 
to infer a more central interaction [^. For events with no dark tracks, no shadow- 
ing is observed while for events with dark tracks, shadowing is enhanced over other 
experimental measurements insensitive to the spatial location. 

The spatial dependence of nuclear shadowing has been studied in charm and 
bottom production While the effect can be important in peripheral collisions 

since nucleons near the nuclear surface are expected to behave as free particles, in 
central collisions where modifications are larger than average, this increase is ~ 1% 
over the average, negligible for our purposes. 



Rates 

All the results presented here are in units of cross section per nucleon pair. The 
total rates in pA and AA interactions are obtained by extrapolation. To calculate the 
rate of charm and bottom production in these interactions as a function of impact 
parameter, the per nucleon QQ cross sections for pA interactions given in this paper 
should be multiplied by the nuclear profile function T^ib), 

N'pih) = a^?TA{b) . (25) 
Likewise, the A A rates are found by multiplying the A A per nucleon cross section by 




Figure 4: The shadowing parameterizations for A = 200 for (a) uy valence quarks, 
(b) u sea quarks, and (c) gluons. The sohd curves show the ratios at fj, — /lo while 
the dashed curves are at = 10 GeV. 



the nuclear overlap function, T^Aib) = J cP sTA{s)TA{\b — s\), 



Nl^ib) = cri^T^b) • (26) 

Thus, integrating over all impact parameters, the pA cross section should scale as A 
times the pp cross section and the AA cross section by times the pp cross section 
without any nuclear effects. Since the only nuclear effect discussed here that can 
change the total cross section is shadowing, total charm and bottom production should 
have a nearly linear A dependence. The exact A dependence is energy dependent 
because higher energies correspond to lower x and larger shadowing effects. One can 
maximize the A dependence by e.g. studying the rapidity distribution of low px QQ 
pairs. 

We note that in the case of pA interactions, the impact parameter dependence of 
charm and bottom production is the same as that of the geometric cross section. This 
similarity arises because the average number of collisions and the average number 
of participants is the same in pA interactions. In nucleus-nucleus collisions, hard 
production is biased toward central collisions, see Ref. for more details. 



Results 

We will show inclusive charm and bottom and exclusive cc and bb distributions at 
several different energies. First we show fixed target charm production at 158, 450, 
and 800 GeV. These energies are chosen because Pb+Pb interactions are studied 
at the CERN SPS at 158 GeV/nucleon and it would be useful to examine charm 
production in pp and pPh interactions at the same energy. However, since this would 
require a secondary beam, we also show the results at the SPS primary proton beam 
energy of 450 GeV. The Fermilab fixed target program uses an 800 GeV proton 
beam which could also be employed in studies of the charm A dependence. Indeed, 
a great many experiments have studied charm production, both at CERN and at 
Fermilab, but unfortunately the statistics for proton beams have generally been poor 
and the newer results have not been used to determine the A dependence except in 
very limited regions of phase space. Better statistics are available for vr beams but 
the data is typically averaged over all targets to enhance statistics rather than to 
study the A dependence. Since the RHIC and LHC heavy ion programs will include 
pA studies, we also show results at the complementary ion-ion energy, ^/s = 200 
GeV/nucleon and 5.5 TeV/nucleon, respectively. Our AA results are only shown for 
RHIC and LHC energies. 

All the pp calculations are shown with {k"^) = 1 GeV^ and all the pA calculations 
are presented with (k^) = 1.35 GeV^ for charm and 2.57 GeV^ for bottom. The AA 
results are calculated with {k'^) = 1.7 GeV^ for charm and 4.16 GeV^ for bottom. We 
note that the (k'^) values used for bottom are somewhat larger than may be suitable 
for the implementation in the MNR code. The pA and AA calculations are modified 



by spatially homogeneous shadowing with the EKS98 parameterization [^, ^ at 
A = 200. 



We first show that the theoretical K factor, i^th, the ratio of the NLO to LO 
cross sections, both calculated with MRST HO in the MNR code is constant in pp 
and pA interactions at 158 GeV in Fig. ^. There is some variation near the edges 
of phase space where the statistics become poor. However, we have checked that 
neither shadowing nor kx broadening leads to significant changes in i^^th- The bare 



quark f^th has also been calculated at collider energies in Ref. |^ with the MRS D— ' 
distributions and was shown to be constant with quark px and xj? as well as pair 
mass and rapidity. The exact value of Kth depends on the set of parton distributions 
chosen and the interaction energy. 

The exclusive and inclusive cross sections are shown for all energies considered in 
Figs. |6|-p!^. The exclusive results are given for the pair p^, rapidity, and distributions 
while the inclusive results are given as a function of quark p^, xp, and rapidity. Note 
that at collider energies, we do not show the xp results since the distributions are 
very narrow at high energies and the rapidity distributions are thus more useful. The 
pair pt and distributions and the quark pt distribution all show the effect of kx 
broadening in pA and AA over pp interactions. The broadening is most apparent 
at fixed target energies and at low pt for the collider energies. The effect is always 
stronger on the pair quantities. On the other hand, information on nuclear shadowing 
is more likely to be obtained from the rapidity distributions. The most useful variable 
for studying shadowing effects is the pair rapidity distribution since Fig. || showed 
that neither kp broadening nor fragmentation affects the pair rapidity. While the 
quark rapidity distribution in Fig. ^ did show a difference between the bare quark 
rapidity distribution and the distributions with kx broadening included, the change 
is due to fragmentation rather than the kp effect. Since shadowing is not symmetric 
around ?/ = 0, we present results over all rapidity and xp- Tables of the pair and 
single quark pp and pair and single quark rapidity for the pp, pA and AA results are 
also given. 

As the center of mass energy is increased, the effect of kp broadening on the cc 
pair Pt and distributions is seen to diminish in Figs. ^ |], [l^, 12, and 16. At 



the same time, the pA pair rapidity distributions exhibit a depletion at forward xp 
relative to the pp distribution. This effect, most obvious on the linear scale of Fig. |1^, 
is due to shadowing and is biased toward forward rapidities where the smallest target 
X values are probed, X2 ~ 3 x 10~^ for pairs with j9t ~ and y = 5. The bb pair 



results in Figs. |IJ and |18| exhibit the same basic trends. Note however that the 



bb pairs tend to have much larger pp than the charm quarks. In AA collisions, the 
rapidity distributions are again symmetric about y = because now negative rapidity 
corresponds to low xi and the shadowing results are thus essentially the product of 
the target shadowing {X2) and the projectile shadowing (xi). 

The pair pp is generally higher than for single quarks at any given energy. At 
NLO, large pp QQ pairs come from e.g. gg — > gg jets where one final state gluon 
splits into QQ, g — > QQ. In this way, the pair gets the entire pp of the gluon 
while the single quark gets only a fraction of the gluon pp. Conversely the single 
quark rapidity is generally greater than that of the pair. That is because i/gg = 
1/2 ln((E<5 + Eq+ p^q + p^q)/ {Eq + Eq- p^Q - p^g)) so that if p^Q is large and 
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Figure 5: Ratio of NLO to LO c quark production in pp and pA interactions at 158 
GeV as a function oi pr, xp, and y. The crosses are the pp results while the circles 
are the pA results. 



positive, it may be balanced by a large negative p^q. The single quark nuclear effects 
show the same trends as those of the pair. Note that the single quark pt distributions 
are all for xp > ^ only. 

Figure [2^ shows the effect of shadowing alone on the pA/pp ratios of the bare 
charm quark distributions at 158 GeV. Without broadening, the ratio decreases 
with pt since X2 ~ 0.5 at y = 0, well into the EMC region where the shadowing 
ratio decreases with x. The effect at low pt can be observed in the xp and y ratios 
where the ratios tend to be greater than unity since X2 is in the antishadowing region, 
X2 ~ 0.14 at y = and ~ 0.05 at y = 1, still in the antishadowing region for gluons 
with low yU^. The LO and NLO ratios, both calculated with MRST HO, are quite 
similar so that the order of the calculation is not strongly influenced by shadowing. 

Since the differences between the distributions are not always obvious in Figs. ^ 



p!9| , we have also calculated the pA to pp ratios in Figs. pi| - |30| and [35| - |38| and the 
AA to pp ratios at the same energies in Figs. |^T[ - |5^ and The pair and single 

quark p^ ratios reflect the effects of the kp broadening, as do the xp distributions 
for \xf\ > 0.5. The kr broadening follows the energy dependent trend predicted 



in Ref. ||3^. The pA/pp ratio increases with px at low energies and, as the energy 
increases, begins to turn over at pair pt ~ 3.5 — 4 GeV for charm, decreasing to unity 
at large pair px- This turnover is at ~ 5 GeV for bottom. As the energy increases, 
shadowing becomes more important at low p^, decreasing the ratio below unity in this 
region. The magnitude of this decrease is a good indication of the strength of gluon 
shadowing, particularly for the pair. While the single quark ratios follow the trends 
of the pair ratios at low p^, the lower statistics of the Monte Carlo at high pt lead to 
large fluctuations. We have checked that taking the single pt ratio for all xp rather 
than Xi? > 0, as is shown in the figures, neither strongly influences the final results nor 
significantly reduces the fluctuations. These fluctuations are particularly strong in cc 
production at y/s = 5.5 TeV where both the pair and single charm pp distributions 
decrease 4-5 orders of magnitude over the pp interval shown. The ratios decrease as 
(f) goes from to vr because the larger kp kick in pA flattens this distribution relative 
to pp so that pA/pp > 1 at ~ and less than unity at ~ vr, compare Fig. ^ to 



Fig. The general trend is similar for all energies because variables more sensitive 
to shadowing are integrated over. The AA/pp ratios exhibit the same trends except 
that the depletion at low pp is stronger and the enhancement is larger before the 
turnover than in the pA/pp ratios. The position of the maximum enhancement is not 
significantly changed. The stronger depletion at low pp is, in part, due to nuclear 
shadowing while the change in the relative enhancement between pA/pp and AA/pp 
is due to the larger (k^) in AA collisions. 

On the other hand, the pA/pp rapidity ratios follow the gluon ratios shown in 
Fig. ^. The peak at x —>■ 1 from Fermi motion appears at large negative rapidity, 
followed by the dip in the EMC region, a broad peak tracing out the antishadowing 
region, and finally decreasing again as the shadowing region is reached at large, pos- 
itive rapidity. The shape is inverted from that in Fig. |^ because large x corresponds 
to large negative rapidity and low x to large positive rapidity. The strength of the 
shadowing effect increases as energy increases while the antishadowing peak, broad 
at 158 GeV, narrows with increasing energy. At collider energies, the full rapidity 



distribution is not shown in the plots so that the x distribution of Fig. ^ is somewhat 
truncated. Both the single and pair rapidity distributions follow the shape of the 
gluon shadowing ratio, suggesting that either can help elucidate gluon modifications 
in the nucleus without significant distortion from kx broadening which affects the 
xp ratios. The AA/pp rapidity ratios show the shadowing in both the target and 
projectile. Now, in addition to the shadowing function for the target nucleus, there 
is a shadowing function for the projectile nucleus that is essentially that of Fig. ^ 
since here low x corresponds to large negative rapidity. Thus the AA/pp ratios are 
symmetric around ?/ = 0. 

Conclusions 

We have studied the effects of kp broadening and nuclear shadowing on charm 
and bottom production over a wide range of energies. We have found that the QQ 
Pt, the single quark px, and the pair distribution are most sensitive to the effects of 
kx broadening while the pair and single quark rapidities are sensitive to shadowing. 
There is significantly less distortion of the nuclear effects in the pair distributions 
which are the most difficult to measure. The same results may be obtained through 
measurements of the single heavy quarks event though the definitions of x and kp are 
smeared relative to those of the pair. 
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Figure 6: Exclusive NLO cc pair production in pp and pA interactions at 158 GeV as 
a function of pt, t/, and (j). The per nucleon cross section is given for pA interactions. 
The dashed curves are the pp results while the solid curves are the pA results. 
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Figure 7: Inclusive NLO c quark production in pp and pA interactions at 158 GeV 
as a function of pr, Xp, and y. The Pt distributions are integrated over xp > only. 
The per nucleon cross section is given for pA interactions. The dashed curves are the 
pp results while the solid curves are the pA results. 
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Figure 8: Exclusive NLO cc pair production in pp and pA interactions at 450 GeV as 
a function of pt, y, and 4>. The per nucleon cross section is given for pA interactions. 
The dashed curves are the pp results while the solid curves are the pA results. 
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Figure 9: Inclusive NLO c quark production in pp and pA interactions at 450 GeV 
as a function of pr, xp, and y. The Pt distributions are integrated over xp > ^ only. 
The per nucleon cross section is given for pA interactions. The dashed curves are the 
pp results while the solid curves are the pA results. 




Figure 10: Exclusive NLO cc pair production in pp and pA interactions at 800 GeV as 
a function of pr, U-i and 0. The per nucleon cross section is given for pA interactions. 
The dashed curves are the pp results while the solid curves are the pA results. 
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Figure 11: Inclusive NLO c quark production in pp and pA interactions at 800 GeV 
as a function oipx, Xp, and y. The Pt distribTitions are integrated over xp > only. 
The per nucleon cross section is given for pA interactions. The dashed curves are the 
pp results while the solid curves are the pA results. 




Figure 12: Exclusive NLO cc pair production in pp, pA and AA interactions at 
^/s = 200 GcV as a function of pr, y, and 0. The per nucleon cross section is given 
for pA and AA interactions. The dot-dashed curves are the pp results, the dashed 
curves are the pA results, and the solid curves are the AA results. 
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Figure 13: Inclusive NLO c quark production in pp, pA and AA interactions at 
^/s — 200 GeV as a function of pr and y. The pr distributions are integrated over 
xp > only. The per nucleon cross section is given for pA and AA interactions. The 
dot-dashed curves are the pp results, the dashed curves are the pA results, and the 
solid curves are the AA results. 








Figure 14: Exclusive NLO bb pair production in pp, pA and AA interactions at 
^/s = 200 GcV as a function of pr, y, and (f). The per nucleon cross section is given 
for pA and AA interactions. The dot-dashed curves are the pp results, the dashed 
curves are the pA results, and the solid curves are the AA results. 
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Figure 15: Inclusive NLO b quark production in pp, pA and AA interactions at 
^/s — 200 GeV as a function of pr and y. The pr distributions are integrated over 
xp > only. The per nucleon cross section is given for pA and AA interactions. The 
dot-dashed curves are the pp results, the dashed curves are the pA results, and the 
solid curves are the AA results. 
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Figure 16: Exclusive NLO cc pair production in pp, pA and A^l interactions at 
^/s = 5.5 TcV as a function of pr, y, and 0. The per nucleon cross section is given 
for pA and AA interactions. The dot-dashed curves are the pp results, the dashed 
curves are the pA results, and the solid curves are the AA results. 
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Figure 17: Inclusive NLO c quark production in pp, pA and AA interactions at 
y/s — 5.5 TeV as a function of px and y. The px distributions are integrated over 
xp > only. The per nucleon cross section is given for pA and AA interactions. The 
dot-dashed curves are the pp results, the dashed curves are the pA results, and the 
solid curves are the AA results. 
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Figure 18: Exclusive NLO bb pair production in pp, pA and AA interactions at 
^/s = 5.5 TcV as a function of pr, y, and 0. The per nucleon cross section is given 
for pA and AA interactions. The dot-dashed curves are the pp results, the dashed 
curves are the pA results, and the solid curves are the AA results. 
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Figure 19: Inclusive NLO b quark production in pp, pA and AA interactions at 
y/s — 5.5 TeV as a function of px and y. The px distributions are integrated over 
xp > only. The per nucleon cross section is given for pA and AA interactions. The 
dot-dashed curves are the pp results, the dashed curves are the pA results, and the 
solid curves are the AA results. 
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Figure 20: Ratio of c quark production in pA to pp interactions at 158 GeV as a 
function of pt, xp, and y. The crosses are the LO results while the circles are the 
NLO results. 




Figure 21: The ratio of pA to pp exclusive NLO cc pair production at 158 GeV as a 
function of pr, y, and 4>. 
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Figure 22: The ratio of pA to pp inclusive NLO c quark production at 158 GeV as a 
function of pr, xp, and y. 
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Figure 23: The ratio of pA to pp exclusive NLO cc pair production at 450 GeV as a 
function of pr, y, and 4>. 
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Figure 24: The ratio of to pp inclusive NLO c quark production at 450 GeV as a 
function of pr, xp, and y. 
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Figure 25: The ratio of pA to pp exclusive NLO cc pair production at 800 GeV as a 
function of pr, y, and 4>. 
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Figure 26: The ratio of pA to pp inclusive NLO c quark production at 800 GeV as a 
function of pr, xp, and y. 




Figure 27: The ratio of pA to pp exclusive NLO cc pair production at ^/s — 200 GeV 
as a function of pt, y, and 4>. 




Figure 28: The ratio of pA to pp inclusive NLO c quark production at ^/s = 200 GeV 
as a function of pr and y. 




Figure 29: The ratio of pA to pp exclusive NLO bb pair production at a/s — 200 GeV 
as a function of pt, y, and 4>. 
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Figure 30: The ratio of pA to pp inclusive NLO b quark production at ^/s = 200 GeV 
as a function of pr and y. 
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Figure 31: The ratio of AA to pp exclusive NLO cc pair production at ^/s — 200 GeV 
as a function of pt, y, and 4>. 
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Figure 32: The ratio of AA to pp inclusive NLO c quark production at y/s = 200 
GeV as a function of pr and y. 




Figure 33: The ratio of AA to pp exclusive NLO bb pair production at a/s — 200 GeV 
as a function of pt, y, and 4>. 
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Figure 34: The ratio of AA to pp inclusive NLO b quark production at y/s = 200 
GeV as a function of pr and y. 




Figure 35: The ratio of pA to pp exclusive NLO cc pair production at ^/s — 5.5 TeV 
as a function of pt, y, and 4>. 
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Figure 36: The ratio of pA to pp inclusive NLO c quark production at ^/s = 5.5 TeV 
as a function of pr and y. 




Figure 37: The ratio of pA to pp exclusive NLO bb pair production at a/s — 5.5 TeV 
as a function of pt, y, and 4>. 
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Figure 38: The ratio of pA to pp inclusive NLO b quark production at ^/s = 5.5 TeV 
as a function of pr and y. 
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Figure 39: The ratio of AA to pp exclusive NLO cc pair production at ^/s — 5.5 TeV 
as a function of pt, y, and 4>. 
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Figure 40: The ratio of A A to pp inclusive NLO c quark production at ^/s — 5.5 TeV 
as a function of and y. 





dace/ dpT 


(nb/GeV) 


duc/dpT 


nb/GeV) 


Pt (GeV) 


pp 


pA 


pp 


pA 


0.25 


2104.000 


1848.000 


1970.400 


2046.000 


0.75 


3470.000 


3440.000 


1669.400 


1909.400 


1.25 


2048.000 


2504.000 


527.000 


705.200 


1.75 


747.600 


1197.000 


115.180 


185.060 


2.25 


198.080 


421.600 


20.100 


38.720 


2.75 


38.140 


116.940 


3.040 


7.006 


3.25 


5.188 


27.220 


0.391 


1.166 


3.75 


0.788 


3.852 


0.044 


0.135 


4.25 


0.107 


0.867 


0.008 


0.017 



Table 3: The NLO exclusive cc and inclusive single charm pt distributions for a proton 
beam of 158 GeV. The exclusive cc distributions are integrated over all rapidity while 
the inclusive single charm px distributions are integrated over xp > Q only. The 
distributions for pp and pA interactions are both given. The per nucleon cross section 
is given for pA interactions. Recall that the intrinsic (fcf,) is 1 GeV^ for pp interactions 
and 1.35 GeV^ for pA interactions. 




Figure 41: The ratio of AA to pp exclusive NLO bb pair production at a/s — 5.5 TeV 
as a function of pt, y, and 4>. 
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Figure 42: The ratio of AA to pp inclusive NLO b quark production at ^/s = 5.5 TeV 
as a function of pt and y. 





dacc/dpT 


(nb/GeV) 


dajdpr ( 


nb/GeV) 


Pt (GeV) 


pp 


pA 


pp 


pA 


0.25 


8044.000 


7156.000 


7214.000 


7342.000 


0.75 


13546.000 


13544.000 


6704.000 


7416.000 


1.25 


8352.000 


10130.000 


2456.000 


3118.000 


1.75 


3268.000 


5020.000 


642.200 


943.800 


2.25 


939.200 


1862.600 


152.060 


246.200 


2.75 


210.600 


552.400 


29.800 


56.960 


3.25 


40.220 


144.480 


6.476 


11.616 


3.75 


8.492 


27.780 


1.405 


2.552 


4.25 


1.934 


6.662 


0.313 


0.598 


4.75 


0.476 


1.481 


0.078 


0.135 


5.25 


0.128 


0.378 


0.020 


0.030 


5.75 


0.042 


0.078 


0.004 


0.007 


6.25 


0.011 


0.024 


0.001 


0.001 



Table 4: The NLO exclusive cc and inclusive single charm distributions for a proton 
beam of 450 GeV. The exclusive cc distributions are integrated over all rapidity while 
the inclusive single charm p^ distributions are integrated over xp > only. The 
distributions for pp and pA interactions are both given. The per nucleon cross section 
is given for pA interactions. Recall that the intrinsic (k^) is 1 GeV^ for pp interactions 
and 1.35 GeV^ for pA interactions. 





dace/ dpT 


(nb/GeV) 


daJdpT (nb/GeV) 


Pt (GeV) 


pp 


pA 


pp 


pA 


0.25 


14168.000 


12530.000 


12486.000 


12300.000 


0.75 


24160.000 


23880.000 


12100.000 


12952.000 


1.25 


15346.000 


18216.000 


4788.000 


5790.000 


1.75 


6266.000 


9312.000 


1383.600 


1899.800 


2.25 


1902.400 


3558.000 


356.200 


543.000 


2.75 


480.600 


1119.000 


86.240 


140.140 


3.25 


109.220 


311.200 


21.080 


34.680 


3.75 


27.720 


71.280 


5.578 


9.570 


4.25 


7.820 


20.280 


1.661 


2.478 


4.75 


2.456 


5.510 


0.493 


0.698 


5.25 


0.787 


1.591 


0.166 


0.245 


5.75 


0.297 


0.491 


0.048 


0.081 


6.25 


0.111 


0.176 


0.018 


0.021 


6.75 


0.044 


0.068 


0.005 


0.007 


7.25 


0.014 


0.027 


0.002 


0.003 



Table 5: The NLO exclusive cc and inclusive single charm p^ distributions for a proton 
beam of 800 GeV. The exclusive cc distributions are integrated over all rapidity while 
the inclusive single charm pt distributions are integrated over xp > only. The 
distributions for pp and pA interactions are both given. The per nucleon cross section 
is given for pA interactions. Recall that the intrinsic (k^) is 1 GeV^ for pp interactions 
and 1.35 GeV^ for pA interactions. 





dacc/dpT (fjh/GeV) 


dac/dpT (/xb/GeV) 


Pt (^breVj 


pp 


pA 


A A 

AA 


pp 


pA 


A A 

AA 


nop; 


lUo.o4U 


o A r\or\ 

84.080 


68.300 


AO A A A 

98.400 


'~7r' Af\ 

75.740 


77.900 


U. ( 


2U4.8UU 


1 A A r\ 
i m.44U 


1 c A Ar\ 

150.640 


1 1 O P f'C\ 

113.660 


A A OOA 

94.221) 


1 AO 1 /I A 

1U2.14U 


i.ZO 


162.540 


158.960 


153.020 


nr\ -1 r>A 

60.160 


55.860 


GO ^ AC\ 

63.540 


i. / 


on o A f\ 

89.340 


101.580 


1 A£? C\ Ad 

106.940 


O tr y< O A 

25.480 


O C A/^ A 

25.960 


OA OOA 

30.880 


o o 

z.zo 


o8.82U 


50.860 


OO A 


1 A OOA 

10.330 


11 A 1 

11. (Jib 


1 A TO A 

14.720 


Z. / 


1 '7/1 A 

15. ^44 


23.080 


OA TA A 

30.700 


A AOO 

4.082 


/I O A/^ 

4.806 


6.076 


C).ZO 


U.Oc/U 


Q 7Qfi 


IzL '^QS 








O. / U 


2.918 


4.246 


7.614 


0.752 


0.902 


1.164 


4.25 


1.433 


2.014 


3.024 


0.377 


0.454 


0.574 


4.75 


0.765 


1.013 


1.433 


0.185 


0.220 


0.285 


5.25 


0.421 


0.543 


0.761 


0.101 


0.123 


0.152 


5.75 


0.243 


0.310 


0.422 


0.058 


0.066 


0.074 


6.25 


0.145 


0.180 


0.233 


0.034 


0.038 


0.048 


6.75 


0.091 


0.110 


0.143 


0.021 


0.026 


0.030 


7.25 


0.058 


0.067 


0.089 


0.013 


0.015 


0.019 


7.75 


0.038 


0.044 


0.055 


0.008 


0.009 


0.010 


8.25 


0.024 


0.029 


0.035 


0.005 


0.006 


0.007 


8.75 


0.017 


0.019 


0.023 


0.004 


0.004 


0.005 


9.25 


0.011 


0.013 


0.016 


0.002 


0.003 


0.003 


9.75 


0.008 


0.009 


0.011 


0.002 


0.002 


0.002 



Tabic 6: The NLO exclusive cc and inclusive single charm pt distributions for colli- 
sions at i/s = 200 GeV. The exclusive cc distributions are integrated over all rapidity 
while the inclusive single charm px distributions are integrated over xp > ^ only. 
The distributions for pp, pA and AA interactions are all given. The per nucleon cross 
section is given for pA and A A interactions. Recall that the intrinsic is 1 GeV^ 
for pp interactions, 1.35 GeV^ for pA interactions, and 1.7 GeV^ for AA. 





dacc/dpT (/xb/GeV) 


dac/dpT (/xb/GeV) 


Pt (GeV) 




vA 


AA 




vA 


AA 


0.25 


407 800 


301 400 


271 800 


1123 800 


740.400 


679 600 

KJ 1 fj * KJ KJ KJ 


0.75 


1791 000 


1221.600 


924.400 


1604 400 


1078 800 

-L. \J 1 KJ * KJ \J KJ 


1005 800 

-L. KJ KJ ^ • KJ KJ KJ 


1.25 


2732 000 


1921 200 

-1- -1- • \j \j 


1395 200 


1202 200 


875 000 

KJ 1 KJ'KJKJKJ 


817.400 


1.75 


2464 000 


1932 400 


1531 600 


755 600 


577 000 


549 600 

KJ^^ tj • yj 


2.25 


1615 200 


1492 800 


1262 000 


428 800 


350 800 

KJ KJ \J ■ KJ \J \J 


335 600 

KJKJKJ ■ v_/ yj yj 


2.75 


938 400 


939 000 


875 600 


242.200 


200 400 

KJ KJ • -L KJ KJ 


200 400 

KJ KJ • -L KJ KJ 


3.25 


540 200 


570 800 


559 000 

'.y '.y t_y • 


143.700 


124.140 


122.820 


3.75 


324 200 


347 200 


348 200 


85 380 


76.260 


78.460 


4.25 


207 800 


214 000 


217 600 


53 500 

KJKJ 'Kjyjyj 


48.260 


49.200 


4.75 


133 820 


141.280 


141.760 


36.640 


31.680 


32 280 

KJ ^ ■ KJ\J 


5.25 


92.220 


97.380 


96.480 


23.160 


21.700 


21.580 


5.75 


65.940 


67.480 


68.420 


16.696 


15.368 


16.480 


6.25 


47.900 


48.180 


49.540 


11.192 


10.042 


10.974 


6.75 


35 580 


35 880 


35.480 


8.354 


8.282 


8.742 


7.25 


26.420 


26 820 


27.000 


6.568 


5.788 


5.694 


7.75 


20.400 


20 200 


20 860 


4.900 


4.650 


4.902 


8.25 


16.306 


15.728 


15.904 


3.378 


3.446 


3.664 


8.75 


12.646 


12.958 


12.840 


2.886 


2.468 


2.688 


9.25 


9.750 


10.112 


11.420 


2.316 


1.973 


1.844 


9.75 


7.740 


8.150 


7.996 


1.631 


1.812 


1.769 


10.25 


6.162 


6.538 


6.474 


1.267 


1.411 


1.527 


10.75 


5.046 


5.468 


5.238 


1.159 


1.099 


1.279 


11.25 


4.216 


4.304 


4.246 


1.024 


0.968 


0.846 


11.75 


3.512 


3.798 


3.504 


0.660 


0.733 


0.765 


12.25 


2.850 


2.928 


3.242 


0.719 


0.554 


0.730 


12.75 


2.474 


2.540 


2.566 


0.406 


0.580 


0.554 


13.25 


2.016 


2.096 


1.990 


0.386 


0.350 


0.405 


13.75 


1.881 


1.862 


1.804 


0.430 


0.349 


0.360 


14.25 


1.601 


1.522 


1.562 


0.283 


0.257 


0.299 


14.75 


1.280 


1.296 


1.334 


0.322 


0.241 


0.238 


15.25 


1.200 


1.146 


1.184 


0.203 


0.271 


0.195 


15.75 


1.011 


0.982 


1.022 


0.176 


0.178 


0.261 


16.25 


0.799 


0.905 


0.947 


0.223 


0.164 


0.171 


16.75 


0.752 


0.720 


0.873 


0.117 


0.166 


0.174 


17.25 


0.647 


0.718 


0.695 


0.102 


0.111 


0.114 


17.75 


0.576 


0.566 


0.604 


0.128 


0.104 


0.171 



Table 7: The NLO exclusive cc and inclusive single charm pt distributions for colli- 
sions at = 5.5 TeV. The exclusive cc distributions are integrated over all rapidity 
while the inclusive single charm Pt distributions are integrated over > only. 
The distributions for pp, pA and AA interactions are all given. The per nucleon cross 
section is given for pA and AA interactions. Recall that the intrinsic {k^) is 1 GeV^ 
for pp interactions, 1.35 GeV^ for pA interactions, and 1.7 GeV^ for 74^4. 





da^ji/dpT (nb/GeV) 


dab/dpT (nb/GeV) 


Vt (GeV) 


TOT) 


vA 


AA 




vA 


AA 


0.25 


196 260 

t_y W • \J \J 


134.740 


115.740 


60 320 

KJ KJ * KJ ^ KJ 


59 260 

^ t_y • ^ KJ KJ 


64.460 


0.75 


512.400 


375 200 


326 800 

KJ \J t \^ KJ \J 


152 900 

-L. ^ t kJ KJ KJ 


150 680 

KJ KJ * KJ K^ KJ 


165 800 

KJ KJ * K^ KJ KJ 


1.25 


632 600 


516 800 


461 600 


199 820 

J- %J %J » K^ ^kJ 


200 800 

^kJ KJ • KJKJKJ 


227 600 

1 tKJKjKj 


1.75 


583 400 


547.400 


549 200 


218 200 

J_ KJ ■ \J \J 


220 200 

\j ■ \j \j 


249 800 

s-i^I-tJ • kjkjkj 


2.25 


411.800 


493 600 


531 400 

Kj KJ J- ■ -L \J \J 


205 200 

\J KJ ■ KJ 


213 800 

_1_ KJ ■ KJ 


245 600 

-L \-/ ■ \J KJ KJ 


2.75 


271.200 


406 000 


484.400 


176 080 

-L 1 KJ * KJ Ky KJ 


188 020 

_L KJ * KJ ^ KJ 


224.800 


3.25 


167 880 


294.400 


399 200 

KJ t_y t_y • ^ KJ KJ 


143.840 


158.220 


186 580 

_L KJ KJ * KJ KJ KJ 


3.75 


108 580 


207 400 


311 200 

KJ -1- -1- • ^ KJ \J 


112.120 


125 600 

J-^KJ • KJkJKJ 


152 840 

_J_ KJ • KJ^IkJ 


4.25 


72.620 


138.440 


223 800 

^ ^kj • Kjyjyj 


83.480 


96 660 

tj\j • \j\j\j 


116.940 


4.75 


51.360 


92 200 


164.180 


64.340 


73.340 


90 860 


5.25 


36.100 


62 080 


112.040 


44.880 


53.240 


68 200 

KJ KJ • KJ KJ 


5.75 


26 980 


41.860 


75.280 


32.120 


40.500 


48.460 


6.25 


20 000 


30.180 


52 800 

KJ ^ • KJkJ \J 


24.240 


27.040 


36 680 

KJKJ • KJ KJkJ 


6.75 


15.304 


21.440 


36 080 

KJ \J ■ \J KJ KJ 


16.774 


20.480 


25 580 

KJ ■ KJ KJKJ 


7.25 


11.450 


15.744 


25 620 


11.866 


14.502 


18.294 


7.75 


8.794 


11.774 


17.932 


8.878 


10.784 


13.860 


8.25 


6.684 


8.946 


13.118 


6.260 


7.496 


9.694 


8.75 


5.188 


6.792 


9.788 


4.486 


5.542 


6.958 


9.25 


4.046 


5.260 


7.328 


3.276 


4.020 


4.902 


9.75 


3.178 


4.036 


5.552 


2.472 


2.798 


3.564 


10.25 


2.464 


3.124 


4.256 


1.741 


2.158 


2.460 


10.75 


1.970 


2.408 


3.254 


1.340 


1.575 


1.866 


11.25 


1.572 


1.906 


2.488 


0.891 


1.168 


1.434 


11.75 


1.220 


1.512 


1.939 


0.810 


0.820 


1.022 


12.25 


0.969 


1.209 


1.554 


0.567 


0.691 


0.835 


12.75 


0.772 


0.929 


1.229 


0.448 


0.503 


0.577 


13.25 


0.625 


0.759 


0.963 


0.301 


0.432 


0.430 


13.75 


0.511 


0.601 


0.756 


0.276 


0.263 


0.347 


14.25 


0.413 


0.478 


0.596 


0.187 


0.210 


0.277 


14.75 


0.342 


0.388 


0.475 


0.166 


0.185 


0.182 


15.25 


0.266 


0.318 


0.389 


0.123 


0.152 


0.162 


15.75 


0.217 


0.250 


0.311 


0.112 


0.088 


0.118 


16.25 


0.176 


0.201 


0.250 


0.043 


0.079 


0.081 


16.75 


0.142 


0.164 


0.195 


0.060 


0.070 


0.077 


17.25 


0.116 


0.134 


0.160 


0.051 


0.051 


0.059 


17.75 


0.095 


0.108 


0.128 


0.033 


0.037 


0.042 



Table 8: The NLO exclusive bb and inclusive single bottom px distributions for colli- 
sions at ^/s = 200 GeV. The exclusive bb distributions are integrated over ah rapidity 
while the inclusive single bottom pt distributions are integrated over xp > only. 
The distributions for pp, pA and AA interactions are all given. The per nucleon cross 
section is given for pA and A A interactions. Recall that the intrinsic (k^) is 1 GeV^ 
for pp interactions, 2.57 GeV^ for pA interactions, and 4.16 GeV^ for AA. 
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Q35 


QQ4 


344 


307 


313 


15.25 


0.775 


0.831 


0.883 


0.302 


0.307 


0.345 


15.75 


0.698 


0.741 


0.778 


0.261 


0.270 


0.276 


16.25 


0.645 


0.680 


0.720 


0.202 


0.232 


0.237 


16.75 


0.572 


0.605 


0.635 


0.207 


0.180 


0.211 


17.25 


0.518 


0.550 


0.581 


0.183 


0.182 


0.205 


17.75 


0.464 


0.485 


0.510 


0.156 


0.166 


0.137 



Table 9: The NLO exclusive bb and inclusive single bottom distributions for colli- 
sions at y/s = 5.5 TeV. The exclusive bb distributions are integrated over all rapidity 
while the inclusive single bottom pr distributions are integrated over xp > only. 
The distributions for pp, pA and AA interactions are all given. The per nucleon cross 
section is given for pA and AA interactions. Recall that the intrinsic is 1 GeV^ 
for PP interactions, 2.57 GeV^ for pA interactions, and 4.16 GeV^ for AA. 





dacc/dy (/xb) 


dac/dy (/ib) 


y 


pp 


pA 


pp 


pA 


-i.ob 


0.000 


n nnn 

0.000 


0.017 


0.014 


-i.oo 


A AAA 

0.009 


A A A O 

0.008 


A AO A 

0.084 


A A '70 

0.072 


-1.45 


A A/^ A 

0.069 


A A/^ A 

0.060 


A C\ A'T 

0.2A7 


0.223 


1 of; 


A O /I O 

0.24d 


A 01 O 

0.218 


A O A 

0.530 


A [X A 

0.505 


1 of; 
-i.UO 


A A 1 

0.601 


A r A 

0.569 


A O O O 

0.888 


A AAA 

0.900 


n o f; 


1 . lOl 


1 1 '^9 

1 . lOZ 


i .zyo 


i .o^y 


-f) fi5 


1.707 


1.781 


1.699 


1.841 


-0.45 


2.237 


2.431 


2.009 


2.230 


-0.25 


2.642 


2.951 


2.254 


2.541 


-0.05 


2.834 


3.231 


2.346 


2.674 


0.05 


2.828 


3.252 


2.340 


2.682 


0.25 


2.644 


3.063 


2.248 


2.593 


0.45 


2.238 


2.602 


2.011 


2.309 


0.65 


1.710 


1.968 


1.692 


1.936 


0.85 


1.129 


1.282 


1.300 


1.462 


1.05 


0.600 


0.655 


0.885 


0.988 


1.25 


0.243 


0.248 


0.532 


0.584 


1.45 


0.069 


0.066 


0.250 


0.261 


1.65 


0.009 


0.008 


0.083 


0.086 


1.85 


0.000 


0.000 


0.017 


0.017 



Table 10: The NLO exclusive cc and inclusive single charm y distributions for a proton 
beam of 158 GeV. The distributions for pp and pA interactions are both given. The 
per nucleon cross section is given for pA interactions. Recall that the intrinsic (k^) 
is 1 GeV^ for pp interactions and 1.35 GeV^ for pA interactions. 





dacc/dy (fih) 


dac/dy {fih) 


y 


VP 


pA 


pp 


„ A 

pA 


-Z.Zo 


f\ f\f\ A 

0.004 


A AA A 

0.004 


A A o r 

0.085 


0.066 


-z.Uo 


0.04d 


A AO O 

0.038 


A OA"! 

0.301 


A O/^O 

0.262 


-i.oo 


0.210 


A 1 AO 

0.192 


A 'VAO 

0.708 


A /I O 

0.648 


1 

-i.DO 


0.604 


0.566 


1.365 


1.346 


-i.4o 


1.367 


1.326 


2.214 


O O o 

2.243 


-i.ZO 


O O O A 

2.389 


2.437 


o 1 r" o 

3.158 


O O A 

3.360 


-i.Uo 


o '^7-^ A 

3.719 


O A/^ A 

3.960 


4.162 


A r 'vo 

4.572 


-U.OO 


5.147 


5.639 


r OA A 

5.200 


r o 'vA 

5.870 


-U.DO 


6.516 


7.329 


6.102 


7.002 


-U.4o 


7.712 


8.813 


6.783 


i-T OO 

7.824 


-U.zo 


8.536 


A AO A 

9.930 


O t A 

7.319 


8.505 


-U.Uo 


8 QQ8 

o.yyo 


iU.OOU 


( .OUD 


O. ( ou 




8.998 


10.540 


7.529 


8.736 


0.25 


8.534 


9.964 


7.284 


8.406 


0.45 


7.702 


8.851 


6.807 


7.761 


0.65 


6.507 


7.344 


6.071 


6.840 


0.85 


5.152 


5.632 


5.193 


5.743 


1.05 


3.722 


3.925 


4.178 


4.543 


1.25 


2.387 


2.402 


3.149 


3.327 


1.45 


1.367 


1.295 


2.176 


2.254 


1.65 


0.605 


0.548 


1.394 


1.405 


1.85 


0.211 


0.183 


0.710 


0.686 


2.05 


0.043 


0.035 


0.292 


0.272 


2.25 


0.004 


0.003 


0.083 


0.077 



Table 11: The NLO exclusive cc and inclusive single charm y distributions for a proton 
beam of 450 GeV. The distributions for pp and pA interactions are both given. The 
per nucleon cross section is given for pA interactions. Recall that the intrinsic {k^) 
is 1 GeV^ for pp interactions and 1.35 GeV^ for pA interactions. 





dacc/dy (fih) 


dajdy (/xb) 


y 


VP 


pA 


VP 


pA 


-z.4o 


0.015 


A A1 O 

0.013 


A 1 OT 

0.187 


0.154 


-Z.Zd 


0.106 


A A A 1 

0.091 


A v r" v 

0.555 


A r" 1 A 

0.510 


-z.Ud 


r\ o r\ A 

0.d90 


A O /I '7 

0.347 


1 1 AA 

1.199 


1 1 AA 

1.109 


-i.oo 


1.019 


A A /I 1 

0.941 


2.164 


o 1 r r 

2.155 


-i.DO 


2.093 


2.042 


o o r» A 

3.320 


O A AA 

3.490 


-i.4o 


3.506 


3.576 


4.685 


r" AA/^ 

5.096 


-i./O 


5.345 


5.741 


1 OA 

6.120 


6.789 


-i.Uo 


7.181 


7.967 


'~7 O A 

7.684 


8.666 


-U.oo 


9.321 


1 A r O A 

10.580 


A A O 1 

9.081 


1 A OAA 

10.390 


-U.DO 


ll.ObO 


1 O O O A 

12.880 


10.340 


1 O AOA 

12.030 


-U.4o 


12.860 


1 A r\t A 

14.910 


1 1 OAA 

11.290 


1 O 1 1 A 

13.110 


nop; 


1/1 A r" A 

14.050 


1 /I OA 

16.420 


11 AO A 

11.980 


1 O A r" A 

13.950 


-U.Ud 


1/1 1 A 

14.610 


1 A /I A 

17.040 


1 O OAA 

12.300 


1/1 OO A 

14.230 


U.Uo 


1/1 /I A 

14.640 


1 'V AAA 

17.020 


1 O O'VA 

12.270 


1/1 1 1 A 

14.110 


U.ZO 


1 /I n^n 


io.iuu 


ii.y / u 


io.DiU 




12.860 


14.420 


11.330 


12.620 


0.65 


11.060 


12.230 


10.320 


11.340 


0.85 


9.324 


9.872 


9.030 


9.675 


1.05 


7.191 


7.314 


7.720 


8.143 


1.25 


5.347 


5.209 


6.141 


6.242 


1.45 


3.510 


3.196 


4.672 


4.737 


1.65 


2.103 


1.819 


3.297 


3.255 


1.85 


1.014 


0.844 


2.189 


2.118 


2.05 


0.391 


0.313 


1.209 


1.097 


2.25 


0.107 


0.078 


0.550 


0.490 


2.15 


0.015 


0.011 


0.193 


0.161 



Table 12: The NLO exclusive cc and inclusive single charm y distributions for a proton 
beam of 800 GeV. The distributions for pp and pA interactions are both given. The 
per nucleon cross section is given for pA interactions. Recall that the intrinsic (/cf^) 
is 1 GeV^ for pp interactions and 1.35 GeV^ for pA interactions. 





da,r^/dy (/xb) 


da^/dy (fih) 


y 


pp 


pA 


AA 


pp 


pA 


AA 


-4.25 


0.003 


0.003 


0.002 


0.148 


0.109 


0.056 


-3.95 


0.125 


0.105 


0.070 


0.939 


0.772 


0.588 


-3.65 


1.022 


0.925 


0.604 


2.988 


3.223 


2.422 


-3.35 


3.398 


3.263 


2.367 


7.240 


7.515 


5.610 


-3.05 


8.264 


8.325 


5.900 


13.420 


15.310 


10.730 


-2.75 


15.050 


16.570 


11.480 


20.550 


22.350 


17.740 


-2.45 


23.850 


26.750 


19.320 


28.580 


33.120 


24.520 


-2.15 


33.690 


38.640 


29.170 


36.600 


43.300 


34.830 


-1.85 


43.780 


50.660 


41.070 


45.650 


52.460 


44.450 


-1.55 


54.080 


63.880 


52.950 


53.330 


61.470 


51.050 


-1.25 


63.830 


73.310 


60.100 


60.650 


67.530 


58.730 


-0.95 


69.960 


78.450 


69.490 


65.320 


70.720 


63.590 


-0.65 


75.700 


82.120 


75.790 


68.690 


72.230 


68.490 


-0.35 


79.700 


82.990 


79.100 


71.950 


72.470 


70.220 


-0.05 


80.720 


80.510 


79.630 


72.860 


70.470 


69.430 


0.05 


80.840 


78.680 


79.990 


72.730 


70.150 


70.350 


0.35 


79.520 


74.340 


79.610 


72.890 


68.970 


69.670 


0.65 


75.930 


67.520 


75.420 


69.340 


61.990 


67.540 


0.95 


70.440 


60.230 


69.290 


65.230 


57.870 


63.320 


1.25 


64.000 


52.670 


59.680 


61.030 


51.720 


59.100 


1.55 


54.240 


43.760 


53.000 


53.110 


44.530 


52.360 


1.85 


43.610 


33.690 


40.990 


46.350 


37.890 


43.830 


2.15 


33.430 


24.880 


29.270 


37.400 


29.630 


35.640 


2.45 


23.680 


17.420 


19.480 


27.760 


21.260 


25.230 


2.75 


15.050 


10.820 


11.600 


20.230 


15.600 


17.760 


3.05 


8.150 


5.742 


5.894 


12.770 


9.649 


10.740 


3.35 


3.403 


2.252 


2.345 


7.394 


5.295 


5.571 


3.65 


1.006 


0.674 


0.606 


2.844 


2.169 


2.426 


3.95 


0.125 


0.090 


0.070 


0.980 


0.662 


0.537 


4.25 


0.003 


0.002 


0.002 


0.180 


0.099 


0.050 



Table 13: The NLO exclusive cc and inclusive single charm y distributions for col- 
hsions at ^/s = 200 GeV. The distributions for pp, pA and AA interactions are all 
given. The per nucleon cross section is given for pA and AA interactions. Recall that 
the intrinsic {k^) is 1 GeV^ for pp interactions, 1.35 GeV^ for pA interactions, and 
1.7 GeV^ for AA. 





da 


cc/dy [jih) 


d(Tc/dy (/ib) 


y 


pp 


pA 


AA 


pp 


pA 


AA 


-4.85 


117.0 


128 


4 


190.7 


186.8 


200 


4 


216.2 


-4.55 


234.1 


278 


6 


265.9 


237.5 


268 


7 


228.0 


-4.25 


369.4 


428 


5 


308.5 


331.2 


346 


9 


291.3 


-3.95 


485.5 


495 


9 


362.1 


426.5 


457 


1 


336.4 


-3.65 


532.9 


586 


9 


367.1 


522.0 


522 


1 


343.8 


-3.35 


561.1 


580 





377.3 


542.4 


553 


1 


364.0 


-3.05 


594.1 


586 


8 


399.1 


590.5 


564 


2 


380.4 


-2.75 


627.1 


602 


2 


407.0 


614.5 


606 


1 


386.3 


-2.45 


650.1 


616 


5 


409.4 


650.6 


565 


6 


396.7 


-2.15 


683.6 


600 





422.8 


642.2 


577 


1 


411.1 


-1.85 


685.5 


587 


8 


428.7 


681.9 


582 


8 


415.0 


-1.55 


710.3 


603 


6 


421.4 


694.6 


588 


5 


426.7 


-1.25 


705.9 


618 


8 


424.9 


672.0 


578 


2 


416.3 


-0.95 


711.8 


595 


5 


435.7 


703.2 


589 


7 


414.9 


-0.65 


735.9 


611 


1 


435.5 


736.7 


617 


2 


435.9 


-0.35 


743.6 


590 


1 


428.8 


692.5 


577 


4 


402.3 


-0.05 


720.6 


556 


8 


414.6 


716.7 


560 


5 


410.3 


0.05 


715.6 


548 





419.2 


721.1 


557 


8 


429.4 


0.35 


741.4 


564 





433.3 


699.3 


552 


3 


427.0 


0.65 


744.0 


567 


4 


436.3 


706.2 


555 


3 


440.3 


0.95 


708.6 


541 


4 


430.2 


701.5 


530 


4 


424.2 


1.25 


702.3 


546 


4 


426.7 


704.2 


520 


6 


403.4 


1.55 


713.9 


506 


9 


423.0 


664.2 


497 


5 


419.4 


1.85 


689.1 


482 


1 


429.9 


678.2 


477 


1 


416.1 


2.15 


687.6 


476 


7 


418.2 


639.9 


462 


9 


396.7 


2.45 


651.9 


470 


5 


409.3 


620.5 


447 


6 


393.0 


2.75 


625.8 


437 


6 


407.2 


606.6 


422 


2 


390.0 


3.05 


589.7 


412 


5 


397.3 


591.6 


415 


5 


382.4 


3.35 


564.0 


384 


7 


374.7 


555.1 


388 


1 


380.9 


3.65 


532.9 


379 


3 


362.5 


516.0 


350 


9 


355.4 


3.95 


479.2 


304 





358.4 


430.2 


297 


4 


325.7 


4.25 


371.3 


257 


4 


321.1 


357.9 


226 


2 


294.6 


4.55 


230.3 


153 


8 


263.9 


226.4 


164 


3 


244.7 


4.85 


112.1 


69.4 


192.8 


193.4 


104 


4 


198.5 



Table 14: The NLO exclusive cc and inclusive single charm y distributions for col- 

hsions at ^/s = 5.5 TeV. The distributions for pp, pA and AA interactions are all 
given. The per nucleon cross section is given for pA and AA interactions. Recall that 
the intrinsic (k^) is 1 GeV^ for pp interactions, 1.35 GeV^ for pA interactions, and 
1.7 GeV^ for AA. 





d^hh/dy (nb) 


daiy/dy (nb) 


y 


pp 


pA 


AA 


pp 


pA 


AA 


-2.65 


0.8 


0.7 


0.5 


14.1 


12.1 


12.0 


-2.45 


4.8 


4.2 


3.2 


33.0 


26.6 


28.7 


-2.25 


16.6 


14.3 


12.1 


61.1 


59.2 


59.1 


-2.05 


41.1 


38.3 


32.2 


103.7 


100.0 


103.2 


-1.85 


82.8 


80.2 


70.4 


150.0 


154.8 


158.6 


-1.65 


138.2 


140.8 


132.5 


203.4 


220.3 


232.5 


-1.45 


213.3 


222.4 


214.9 


262.8 


287.3 


302.5 


-1.25 


289.1 


313.4 


321.9 


323.3 


361.1 


399.4 


-1.05 


382.7 


417.3 


441.6 


383.9 


432.0 


485.2 


-0.85 


463.7 


518.7 


576.3 


443.2 


489.1 


558.3 


-0.65 


543.8 


609.8 


702.0 


482.7 


548.9 


627.4 


-0.45 


606.8 


692.5 


814.8 


520.3 


589.4 


681.2 


-0.25 


651.4 


739.2 


890.1 


547.5 


615.2 


722.1 


-0.05 


674.7 


764.2 


930.9 


555.3 


617.2 


738.3 


0.05 


674.1 


762.7 


933.5 


553.5 


615.6 


734.7 


0.25 


650.8 


726.9 


890.2 


545.0 


606.1 


725.8 


0.45 


607.3 


671.0 


814.3 


520.1 


570.5 


684.9 


0.65 


543.3 


582.9 


701.9 


481.7 


522.5 


632.2 


0.85 


463.9 


491.0 


576.4 


440.1 


466.3 


559.4 


1.05 


383.1 


391.0 


442.2 


387.5 


406.7 


481.3 


1.25 


289.0 


291.1 


322.1 


323.5 


338.0 


399.0 


1.45 


213.3 


205.7 


215.1 


264.0 


269.3 


304.8 


1.65 


138.7 


130.4 


132.3 


205.7 


210.2 


233.8 


1.85 


82.9 


75.3 


70.2 


148.6 


147.1 


159.0 


2.05 


41.2 


36.4 


32.1 


102.1 


99.0 


104.3 


2.25 


16.6 


13.7 


12.1 


59.9 


58.3 


58.7 


2.45 


4.8 


4.0 


3.2 


31.3 


30.7 


28.2 


2.65 


0.8 


0.6 


0.5 


13.5 


12.8 


11.4 



Tabic 15: The NLO exclusive hb and inclusive single bottom y distributions for col- 
lisions at ^/s = 200 GeV. The distributions for pp, pA and AA interactions are all 
given. The per nucleon cross section is given for pA and AA interactions. Recall that 
the intrinsic (k^) is 1 GeV^ for pp interactions, 2.57 GeV^ for pA interactions, and 
4.16 GeV^ for AA. 







dcTb/dy (/xb) 


y 


pp 


pA 


AA 


pp 


pA 


AA 


-4.85 


1.811 


1.831 


1.189 


3.012 


3.239 


2.484 


-4.55 


2.940 


3.117 


2.157 


4.652 


5.331 


3.699 


-4.25 


5.035 


5.625 


3.782 


6.328 


7.331 


5.407 


-3.95 


7.070 


7.946 


5.609 


8.606 


9.596 


7.215 


-3.65 


9.872 


11.200 


8.283 


10.500 


12.110 


9.164 


-3.35 


12.260 


13.800 


10.220 


12.780 


14.330 


11.140 


-3.05 


14.330 


16.330 


12.480 


14.770 


16.100 


12.750 


-2.75 


17.350 


19.110 


14.890 


16.950 


18.230 


13.600 


-2.45 


19.100 


21.030 


16.430 


18.730 


20.040 


15.780 


-2.15 


21.350 


22.880 


17.990 


20.070 


20.650 


16.870 


-1.85 


22.890 


23.940 


19.050 


21.300 


22.290 


17.730 


-1.55 


24.010 


23.790 


19.440 


23.250 


22.740 


18.530 


-1.25 


25.320 


24.830 


20.760 


24.040 


23.790 


19.360 


-0.95 


26.500 


25.570 


20.850 


25.360 


23.960 


20.110 


-0.65 


26.660 


25.160 


21.330 


25.380 


24.050 


19.720 


-0.35 


27.460 


25.330 


21.770 


25.480 


23.380 


20.670 


-0.05 


27.050 


24.680 


21.470 


26.170 


23.810 


20.650 


0.05 


27.010 


24.500 


21.500 


25.790 


23.920 


20.390 


0.35 


27.330 


24.130 


21.680 


25.320 


22.810 


20.150 


0.65 


26.660 


23.340 


21.300 


25.380 


21.990 


20.100 


0.95 


26.460 


22.770 


20.900 


25.480 


21.810 


19.870 


1.25 


25.310 


21.310 


20.760 


23.830 


20.230 


19.310 


1.55 


24.080 


19.630 


19.390 


22.540 


19.050 


18.840 


1.85 


22.830 


19.060 


18.930 


21.370 


17.920 


17.640 


2.15 


21.320 


17.380 


18.040 


20.000 


16.690 


16.930 


2.45 


19.100 


15.280 


16.470 


18.770 


15.570 


15.870 


2.75 


17.360 


13.580 


14.940 


17.320 


13.780 


14.260 


3.05 


14.410 


11.260 


12.460 


14.670 


11.580 


12.710 


3.35 


12.140 


9.320 


10.240 


12.740 


10.040 


10.770 


3.65 


9.831 


7.393 


8.249 


10.430 


7.929 


8.860 


3.95 


7.049 


5.294 


5.607 


8.431 


6.380 


7.131 


4.25 


5.017 


3.721 


3.810 


6.190 


4.889 


5.378 


4.55 


2.925 


2.147 


2.171 


4.660 


3.450 


3.950 


4.85 


1.797 


1.290 


1.193 


3.080 


2.222 


2.346 



Table 16: The NLO exclusive bb and inclusive single bottom y distributions for col- 
lisions at a/s = 5.5 TeV. The distributions for pp, pA and AA interactions are all 
given. The per nucleon cross section is given for pA and AA interactions. Recall that 
the intrinsic {k^) is 1 GeV^ for pp interactions, 2.57 GeV^ for pA interactions, and 
4.16 GeV^ for AA. 



